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' Abstract. In this paper, we examine the Shapiro delay caused by the close approach 

■ of an asteroid to the LISA constellation. We find that the probability that such an 

^jrjl event occurs at a detectable level during the time interval of the mission is smaller 

^ ■ than 1 % 
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^ . 1. Introduction 

o\ 

^ ' LISA (Laser Interferometer Space Antenna) [1], a space experiment devoted to the 
detection of gravitational waves, is a nearly equilateral triangular constellation of three 
^ ' spacecraft, the center of mass of which follows an Earth- like orbit. This constellation 
' is located about 20 degrees behind Earth, and the distances between spacecraft are 
planed to be of the order of 5 millions kilometers. The spacecraft exchange optical laser 
beams, and the oscillations of the distances between spacecraft are interferometrically 
monitored. Thanks to this TDI (Time Delay Interferometry) method, gravitational 
waves signatures are tracked in the [~ 10~^, ~ 10~^] Hz frequency domain. To use 
this method, spacecraft inter-distances (and their fluctuations) must be known with a 
precision such that the gravitational field of the Sun has to be modeled in a relativistic 
framework [2]. The analysis of gravitational waves signal provides another useful way 
to explore the universe, giving astrophysical and cosmological information, inaccessible 
from the electromagnetic window and complementary to it. 

The gravitational field of the solar system is generally modeled including the Sun 
and planets only. However, it is known that the terrestrial orbit is frequently crossed 
by asteroids [3,4], refered to as geo-cruisers (GC) in the following. Since LISA is on an 
Earth-like orbit, close encounters with GCs are expected to occur, and the gravitational 
field of a GC passing close to the LISA constellation can generate a signal in the data. 
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In a recent paper, Vinet [5] examined the direct action of a GC on the LISA 
constellation. The author addressed the shift in position of a station due to the direct 
action of the asteroid's gravitational field. He finds that this effect leads to a measurable 
signal if the involved GC passes sufficiently close to the spacecraft concerned. 

In this paper, we are interested in another aspect of the interaction between 
asteroids and LISA. If a GC passes close to the segment joining two spacecraft, its 
gravitational field affects the hght-distance between them by Shapiro delay on the laser 
beam used for that very measurement. According to the close encounter parameters, 
this signal can have a duration such that it falls in the LISA frequency domain. Hence, 
this signal has to be distinguished from the expected extra-solar system gravitational 
wave signal to be tracked in the data. The aim of the present study is to examine if 
such an (impulsive) event is likely to occur. We find that, while the effect is effectively 
measurable for sufficiently massive GCs passing close to the light beam, the probability 
of the occurence of such an event at a detectable level, during the time interval of the 
mission, is quite negligible. 

2. Conditions for an asteroid encounter to cause a relevant Shapiro delay 

Let us consider an GC passing close to the light beam linking two LISA spacecraft A 
and B. The spacetime geometry in which the beam propagates can be formally written 
as 



gravitational field due to the Sun and planets, '' the part due to the close asteroid, 
and h!"^p^'^ the part of the gravitational wave. The gravitational wave term induces a 
change in the distance between the two stations of the order of 



where h is the characteristic amphtude of /i)^,^' ' , and L the distance between A and B. 
On the other hand, h^^^'^ is of the order of 2Gm/(rc^), where m is the asteroid mass. 
The close approach induces a Shapiro time delay 5t in the flight time of the photon, 
hence a change in the light-distance given in [6] 



where ta (resp. tb) is the distance between the GC and spacecraft A (resp. B) and A 
the distance between the GC and the segment joining the spacecraft A and B. Let h be 
the impact parameter of the encounter (minimum value of A during the approach). In 
the case where h « L (we will find that it is a necessary condition for the signal to be 
observable), it is easy to see that 5L(^^*), the maximum possible value for 51, satisfies 



where 77^/3 





^^(ast.) < 
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2.1. Conditions on the amplitude 

Let -ffmin be the smallest value of h accessible to the experiment. The necessary condition 
for the GC to generate a gravitational signal with sufficient amplitude to be detectable 
(see the next sub-section for the necessary condition related to the frequency domain) 
writes 

— — In - > H^^^L. 
& b 

Let p and D be the density and the (mean) diameter of the considered asteroid 
respectively. The above unequality leads to the following condition on the impact 
parameter 

for the asteroid signal to be observable. This gives, numerically, 

b < (5.10^ km) exp | -8000.^^ ( ^) ' ( ' \ (2) 

I 10-20 g/cm3y \lkinj J 

where we have taken L = 5.10^ km, the average inter-distance between spacecraft. Let 
us consider an GC of 10 km in diameter (resp. 15 km). One finds (with p = 2 g/cm^ 
and ifmin = 10~^^) b < 1700 km (resp. 470000 km). We note that for a 8 km diameter 
asteroid, the impact parameter should be smaller than 1 km, that is smaller than the 
asteroid radius, which means that the beam would be occulted. 



2.2. Conditions for the signal to fall in LISA's frequency interval 

The characteristic time of the encounter is given by r ~ b/V, where V is the relative 
velocity of the GC with respect to LISA's center of mass. This means that the 
fundamental frequency, in the Fourier representation of the signal, is of the order of 
V/b. Then, a second necessary condition is that V/b should be inside the frequency 
interval of LISA for the signal to be detectable. Simulations with fictitious impulsive 
signals (i.e. signals with limited duration), the duration of which range from 1 to 10^ s, 
confirm that -ffmin is always > 10-^°^ and that the detection is not efficient outside the 
interval [10 s, 10^ s] (see appendix I). 

Since one should have r in the time interval [10 s, 10^ s], b has to satisfy the 
additional condition 

{lOs).V<b<{10's).V (3) 

besides condition ([2]). Since \^ ~ 15 km/s, ([3]) leads to 

150 km<b< 150 000 km, (4) 

one sees that only a small number of GCs will effectively be relevant to LISA at 
-f^min = 10"2°, verifying both conditions ([2]) and ([3]). Indeed, only asteroids larger 
than 9 km in diameter can generate a signal at a detectable level with characteristic 
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encounter times in this interval. From astronomical observations, only about 15 GCs 
are larger than 9 km in diameter [4]. If one takes p = 2.7 g/cm^, this limit in diameter 
becomes 8 km, and about 20 GCs are larger than 8 km in diameter. 



3. Probability of a relevant encounter 

Let n (> Dq) be the mean number density of GCs with a diameter D > Dq in the 
neighbourhood of the Earth orbit. Let V be the mean relative velocity of GCs and the 
Earth. The number of GCs, of diameter larger than Dq, passing at a distance between 
b and b + db from the segment [A,B] (with b « L) during a time interval dt, is of the 
order of 2n (> Dq) .L.db.V.dt. Let Tlisa be the duration of the LISA mission. From eq. 
([1]), the condition of detectability by LISA is D > Dq, with 

3 _ 3 HL 
° ~ 4^Gpln(L/6)' 

The number E of events observed during the duration of the mission is then of the order 
of 

E ~ GLFTlisa / n{D>Do)db 



mm 



since there are three arms in the LISA configuration. The lower and upper bounds 
6min and 6max are the minimal and maximal values of b, related to the LISA frequency 
sensitivity curve for impulsive events (jlj). The total number of GCs of diameter > Dq 
is estimated to be [4] 



A.(>C„). 1090 (^j . (5) 

To evaluate E, only the density number in the vicinity of Earth is needed, not N, the 
total number of GCs. From and the estimate obtained in the appendix II, the mean 
number of GCs per unit volume (per (AU)^) in the vicinity of the Earth orbit is 

'1.95 



n(> Dq) ~ 94 



Dn 



1 km 



Then 



E 



1 A.U.] 



_^ n 0.65 



In 



db. 



2 g/cm^ V10-2° / b 

We have taken = 15 km/s and Tlisa = 3 yrs. The values of 6min and 6max were given 
in (jlj), and ln(L/6) in the integral varies in the interval [~ 3.5; ~ 10.4]. The minimal 
amplitude -ffmin depends on the characteristic time of the encounter, hence on b, but, 
as stated before, it can be bounded by 10~^°. Since p is always of the order of 2 g/cm^ 
(for asteroids, it belongs to the interval [1.3 g/cm^; 2.7 g/cm^]), the number of relevant 
events during the LISA mission is bounded by 

E < 1.65 10"^ 
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which means that the probabihty to observe one event caused by a Shapiro delay related 
to a close GC approach is quite negligible. This number (probability) becomes 2 10~^ 
if one takes p = 2.7 g/cm?. 

4. Discussion 

The present study leads to the conclusion that GCs will not perturb significantly the 
LISA mission through related Shapiro effect. While it could appear that this (not very 
exciting) result is not a surprise, a dedicated study was required. Indeed, we have shown 
that the close approach of one of the largest geo-cruisers is susceptible to result in a 
detectable signal. The low number of relevant encounters comes from the statistical 
aspect of the problem, but not directly from the physical properties (masses) of asteroids 
neither from the geometry of the possible encounters. 

It is worth pointing out that, since only objects with a diameter larger than 
8 km are relevant, all the corresponding geo-cruisers are known. Hence, if a relevant 
close encounter with LISA occurs, the corresponding geo-cruiser motion will have been 
accurately monitored, in such a way that it should be easy to veto the resulting signal. 
Consequently, it would be useless to make extensive templates of asteroid's signals, in 
order to track such events during the whole LISA mission. 

In reference [5], the direct effect of an asteroid on the motion of a LISA spacecraft 
has been addressed. A statistical analysis of this direct effect would be of interest and 
remains to be made. To achieve this study, an analytical expression relating the impact 
parameter of the encounter and the minimal value of the asteroid's diameter for which 
the direct effect is detectable is required. Such an expression is not exphcitly provided 
in [5], but it appears from the curves on fig. 3 of [5] that asteroids with a size of (say) 
100 m lead to a detectable effect for impact parameters in the interval (HI). This shows 
that the direct effect examined in [5] is considerably more important compared to the 
Shapiro effect examined in the current paper, as the latter requires larger GCs (at least 
8 km in diameter). 

5. Conclusion 

The present study shows that the possibility of detecting an asteroid through Shapiro 
delay by the LISA mission : 

- concerns only a small number of geo-cruisers (about twenty at best). An 
occultation of the laser beam occurs before the detection condition is satisfied for geo- 
cruisers with a diameter less than ~ 8 km ; 

- has a very low probability to occur during the time interval of the mission, at best 
of the order of some 10~^. 
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Appendix I : Minimal detectable amplitude 

The signal to noise ratio averaged over all sky directions and polarizations can be 
expressed as [7] 



where Seff (z^) is the effective sensitivity of LISA. In our calculations, we adopt the 
position noise budget for a standard Michelson configuration, including the contribution 
of the galactic binary WD- WD confusion noise [8] . 

Let h{t) = Hf{t) be a gravitational signal of duration T and of amplitude H (the 
function f{t) being of amplitude unity). The corresponding spectral density (u) can 
be expressed as 

SUi^) = H'\fii^,T)\\ 

f being the Fourier transform of /. 

Combining the above equations, one obtains for the minimal detectable amplitude 



where 



oo 

du 



'0 Seff (u) 

Following the convention adopted in the LISA community, we assumed a detectability 
threshold of (S/N)^.^ = 5. 

For instance, let us consider a (fictitious) signal such that f{t) is zero outside 
the interval [0,T], which is unity inside the interval [^1,^2]? and which is linear in the 
intervals [0,ti] and [^2,7"], in such a way that the whole signal is continuous. For this 
signal. Figure [1] exhibits the minimum detectable value -ffmin as a function of the global 
duration T, for ^2 = T/2 and ti = T/2, T/4, T/8 and T/16 respectively. It confirms 
that ifmin is always > 10~^°, and that this conclusion does not dependent drastically 
on the precise signal's profile. It also shows that the detection is not efficient for a 
duration outside the interval [10 s, 10^ s], i.e. for durations outside the LISA frequency 
interval. More precisely, the detection is not efficient for a duration < 10 s. It can be 
efficient for a duration > 10^ s, but only for a highly non-symmetric signal, for which the 
time derivative takes values significantly larger than H/T. However, when a GC passes 
close to LISA, the relative velocity of the encounter is nearly constant, so that the GC 
gravitational field varies in a very regular and quasi-symmetric way. The corresponding 
signal is such that its time derivative is never significantly larger than /imax/T (^max 
being the maximum value of the signal and r its characteristic duration). Hence the 
reasonable assumption that the detection is not efficient when r is outside the interval 
[10 s,10^ s]. 
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Appendix II : From total asteroid distribution to volumic distribution near 
Earth orbit 

The distribution of GCs with respect to orbital elements (a, e, i) is given in reference 
[3]. Since one is interested in an order of magnitude estimate rather than in precise 
results, let us make the following assumptions and simplifications : 

- the diameter distribution of GCs is independent of the orbital-element 
distributions ; 

- the distribution in inclination is limited to the interval [Ojimax], in which it is 
uniform. 

Let us consider an asteroid with orbital elements a and e. The probability that this 
asteroid is at a distance from the Sun in the interval [r, r + dr], at an arbitrary time, 
is given by dP (r, r + dr) = 2dr/ (9 |f|) , where 9 is the period. Then, f is given by the 
energy integral, and one finds 

dr 



dP (r, r + dr) — ^ 

Let p (a, e) be the density distribution in a and e, so that p (a, e) da.de is the probability 
that a GC, arbitrarily chosen in the population, has its semi-major axis and eccentricity 
in the intervals [a,a + da] and [e,e + de], respectively. Using the assumption on 
inclination distribution, one finds that the number density of GCs on an Earth orbit is 
related to the total population N by 

N f p{a, e) da.de 



n 



I 



47r2sini„... ./{.,e} /e2^2 _ (i _ 



where one has replaced r by unity (1 A.U.). In this expression, a is expressed in A.U., 
n in (A.f/.)~^ and {a, e} is the integration domain in the a — e plane. From [3], let 
us consider that the integration domain is bounded by 0.5 < a < 3 and 0.2 < e < 1. 
Besides, for an asteroid to be a GC, one has necessarily a (1 — e) < 1 < a (1 -|- e). Since 
one is only interested in an order of magnitude, let us repace p{a,e) by its mean value 
(p (a, e)) = (1 + In 1.44)"^ ~ 0.7328 in the integral. One finds 

A-K^ sm 

withX = 37r/10 + arcsin(2/3) + 2/31n[(3 + \/5) /2] 2.3138. Prom [3] , a reasonnable 
value for imax is 30 degrees. This leads to 

n{> Dq) ~ 0.086 A^(> Dq) . 

This is the mean density (per [A.U.)^) of GCs with diameter D > Dq, at one 
astronomical unit from the Sun. 




Figure 1. minimum signal's amplitude detectable as a function of its duration, for 
various fictitious signals with profile described in appendix I. 
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